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udss( , & + -meson) state with J p = 1 + is expected to exist in the mass region 



model, a possible resonant state of udss(J = 1 + ) is suggested to appear at 
around 1.4 GeV with the width 0(20 ~ 80) MeV. We propose that the 
meson is a good candidate for the tetraquark search, which would be observed 
in the K + K + ir~ decay channel. 



■ I. INTRODUCTION 



The possibility of multiquark states has been discussed for a long time [1-11]. In partic- 
ular, the possible qqqq states have been suggested in many theoretical efforts to understand 
light scalar mesons (see, for example, Refs. [1,3,4]). The 4g states were proposed in de- 
scriptions of / (600) and / (980), where the strong attraction between (qq)s and (qq)s plays 
an important role [1,3]. Here, (qq)s and (qq)^ denote the color- anti-triplet quark pair and 
the color-triplet anti-quark pair, respectively. On the other hand, the KK molecule states 
were suggested to understand the properties of /o(980) and a (980) [4]. Since the masses 
of negative-parity mesons are expected to be above 1 GeV in a naive interpretation with 
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P-wave qq states, it is considered that these scalar mesons below 1 GeV may be hybrids of 
P-wave qq and compact (gg) 3 (<?<?) 3 with meson-meson tails in the outer region, as argued 
in Ref. [12]. Even if the 4q components are dominant in a certain meson whose minimal 
content is 2q, it is difficult to find direct evidence of the 4q components due to mixing with 
the conventional qq state via the annihilation of qq pairs. Our main interest here is the 
possibility of narrow "tetraquark" states, whose minimal quark content is 4g. 

The recent observation of P s j(2317) [13] and reports of the pentaquark baryon 
Q + {uudds) [14-22] revived motivation for experimental and theoretical studies on multi- 
quarks in hadron physics, though the existence of 6 + is yet to be well established. One of the 
striking characteristics of the + is its narrow width. For a theoretical interpretation of why 
the + is extremely narrow, the possibility of the spin-parity J p = l/2 + and J p = 3/2~ has 
been discussed by many groups [11,23-29]. Since only the P-wave and P-wave are allowed 
in NK decays from the J p = l/2 + and J p = 3/2~ states, respectively, the width should be 
suppressed due to a high centrifugal barrier. The transition into meson-baryon states should 
be further suppressed if the pentaquark has an exotic color configuration {qq)^{qq)^q- The 
factor 1/3 in the transition appears from the overlap of the color wave functions of quarks. 
In addition to suppression due to the color degrees of freedom of 5 quarks, another sup- 
pression effect can be considered in flux-tube pictures [6,26,30-32], because the transitions 
between different flux-tube topologies are suppressed due to a rearrangement of the gluon 
field. This means that the decays from such exotic flux-tube configurations as {qq)i{qq)iq 
(Fig. 1(e)) into meson-baryon-like (qqq)i(qq)i (Fig. 1(f)) are suppressed. In general, the 
different flux-tube configuration appears in multiquarks that contain more than 3 quarks, as 
shown in Fig. 1, and the coupling between the disconnected tube and the connected tube 
topologies should be strongly suppressed. 

The predicted spin-parity J p = l/2 + and J p = 3/2~ of the pentaquark uudds are abnor- 
mal in a naive quark picture, where the J p = l/2~ should be the lowest state, while other 
spin-parity states are expected to be highly excited. Originally, a narrow J p = l/2 + was 
predicted in the Skyrme soliton model by Diakonov et al. [33]. For a theoretical explanation 
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of the J p = l/2 + state from the point of view of the constituent quarks, the diquark picture 
and the triquark picture are proposed in Refs. [23,24]. Recently, constituent quark model 
calculations [26,27] suggested an abnormal level structure in the compact uudds system with 
the {qq)i{qq)zQ configuration, where the masses of the J p = l/2 + and/or J p = 3/2~ states 
may degenerate with the J p = l/2~ state. 

We now turn to a discussion on the possibility of tetraquarks. By replacing a wd-diquark 
in the G + with an s quark, it is natural to expect that a tetraquark with the udss content 
may exist at nearly the same energy region. In order to search for a narrow tetraquark, we 
follow the analogy of the theoretical explanation why the pentaquark 6 + can be narrow. 
Firstly, one should consider those states with unnatural spin and parity, which cannot decay 
into two light hadrons (psuedscalar mesons) in the S-wave channel. Second, the exotic flux- 
tube configurations with connected tubes (Fig. 1(c) and (e)) would be essential to stabilize 
the exotic hadrons. We then propose a J p = 1 + udss state with the (qq) 3(^)3 configuration 
as a candidate of narrow tetraquark states. It should be stressed that two-body KK decays 
from any J p = 1 + udss state are forbidden because of conservation of the total spin and 
parity. The lowest threshold energy of the allowed two-body decays is 1.39 GeV for the 
KK* (895) channel. If the mass of the J p = 1 + udss state lies below the KK*, two-meson 
decay channels are closed, and hence its width must be narrow. 

The 4q states have been proposed by Jaffe in 1977 [1]. The (gg)g diquark and the (qq)s 
anti-diquark correlations play an important role in the stability of the 4g state, because 
there exists attraction between (qq)^ and (qq)^ due to the confining and the one-gluon- 
exchange(OGE) potential. Moreover, it is known that the spin-zero flavor- singlet (qq)s 
diquark is favored, because it gains the color-magnetic interaction. Thus, the diquark(anti- 
diquark) correlation leads to confined 4q states with the {qq)i{qq)z color configuration, which 
might make the Aq system compact and stable. In the J p = 1 + state of the tetraquark 
{ud)z(ss)z, one can naturally expect that the lowest is the isoscalar J p = 1 + state with a 
spin-zero (ud)^ and a spin-one (ss)^ in a spatially symmetric orbit. In the spatially symmetric 
(ss)s, the spin-zero(spin-singlet) configuration is forbidden, and hence the (ss)s must have 
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spin-one. Although the spin-one (ss)s feels some repulsive color-magnetic interaction, the 
repulsion is expected to be small because the color-magnetic term in the OGE potential is 
suppressed by the quark-mass factor m~ 2 . In the flux-tube model, the (qq) 3(^)3 state has 
the exotic tube topology shown in Fig. 1(c). Therefore, its coupling with two-hadron states 
should be small, due to the suppressed transitions between the different tube topologies, (c) 
and (d). 

The tetraquark udss states are discussed in Ref. [1], and noted as E( KK ymewns. In 
Ref. [1], the theoretical mass for the isoscalar udss(J p = 1 + ) state is predicted to be 1.65 
GeV in the MIT bag model [1]. Recently, the isoscalar udss in the flavor 10 group was 
suggested in analogy with the + by T. Burns et al. [34] and by Karliner and Lipkin [47]. 
In Ref. [47], the possibility of a + state is discussed. However, we should remark that 
the J p = + is not allowed in the isoscalor udss system within the spatially symmetric 
configuration, and hence, the J p = + is expected to be unfavored. The tetraquark udss is 
called a -# + -meson in Ref. [34], where the J p = 1~ state with the orbital angular momentum 
L — 1 is predicted in the mass region ~ 1.6 GeV. Although the J p = 1~ state may gain 
color-magnetic attraction, it needs the L = 1 excitation energy, and is expected to be higher 
than the spatially symmetric state. Another claim for the -# + (J p = 1 _ ) state is that it can 
decay into P-wave KK states. The centrifugal barrier may not be high enough to stabilize 
the state much above the threshold energy. Therefore, we think that the f) + (J p = 1 + ) is a 
better candidate for narrow tetraquarks. 

In this paper, we consider the -$ + -meson with J p = 1 + by a constituent quark model. The 
theoretical method of the calculations is the same as that applied to the pentaquark study 
in Ref. [26]. Namely, we apply the flux-tube quark model with antisymmetrized molecular 
dynamics(AMD) [35,36] to the 4g systems. Based on the picture of a flux-tube model, we 
ignore the coupling between configurations shown in Fig. 1(c) and Fig. 1(d), and solve the 
4g dynamics with the variational method in the model space {qq)^{qq)^ shown in Fig. 1(c). 
The Coulomb and color-magnetic terms of the OGE potential and the string potential are 
taken into account. In order to evaluate the $ + (J p = 1 + ) mass, we adopt the observed 6 + 
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mass as an input as well as the normal hadron spectra. We also try to interpret a /i-meson 
in 1.4 ~ 1.6 GeV region with the 4 q state, which would help to check the reliability of the 
present calculations. The widths of these states are also discussed. 



II. HAMILTON! AN 

The adopted Hamiltonian is the same as that of previous work [26] as H — H + Hi + 
Hf, where H consists of the mass and kinetic terms, H I represents the short-range OGE 
interaction, and Hf is the string potential given by the energy of the flux tubes. The quarks 
are treated as non-relativistic spin-| Fermions. The OGE potential consists of Coulomb and 
the color-magnetic interactions, as 

1 2tx 

H I = a c Y / FiFj[— ~ oZTZT 8 ^)^ ' (!) 

Here, a c is the quark-gluon coupling constant, and FiFj is defined by Y^a=i,-,sFi t F", where 
is the generator of color 577(3), |Af for quarks and — |(Af)* for anti-quarks, m; is the 
quark mass m q for u and d quarks, and m s for a s quark. The usual Sfcj) function in the 
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spin-spin interaction is replaced by a finite-range Gaussian, s(r,ij) = 2 J^ A exp 

In the flux-tube quark model [6,37], the confining string potential is written as Hf = 
oLf — M°, where a is the string tension, Lf is the minimum length of the flux tubes, and 
M° is the zero-point string energy. For the meson and 3g-baryon systems, the flux tube 
configurations are the linear line and the F-type configuration with a junction, as shown in 
Figs. 1(a) and (b), respectively. For the q 2 <f mesons and q A q baryons, the exotic topologies 
Figs. 1(c) and (e) appear corresponding to the {qq)i{qq)^ an d (<?<?) 3 (<?<?) 3<L m addition to 
the normal two-hadron configurations (Fig. 1(d) and (f)). In principle, besides these color 
configurations, other color configurations are possible in totally color-singlet q 2 q 2 and q A q 
systems by incorporating a color-symmetric (qq)^ pair. However, since such a string from 
the (qq)e is an excited one, and is unfavored in the strong-coupling limit of lattice QCD [38], 
we should consider only color-3 flux tubes as the elementary tubes. The string potentials 
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given by the tube lengths of the configuration Fig. 1(b), (c) and (e) are supported by lattice 
QCD calculations [39-41]. 

In the practical calculation of the expectation values of the string potential (<E»|if/|<E>) with 
respect to a meson state^,^), a three-quark state($g3), a (gg) 3 (<?<?) 3 st&te(§( qq ) 5 (qq) 3 ) , and a 
{qq)^{qq)^,q state ($(99)3(55)59), the minimum length of the flux tubes Lj is approximated by 
a linear combination of two-body distances r\j as, 

L f « rn in (2) 

~ ^(n2 + r 23 +r 3 i) in (^if/l^s), (3) 

L / ~ ^( r i2 + m) + \{rn + rig + r 2l + r 22 ) in ($( w ) 5 (ot) s I#/I$(ot) 3 (m)3)> (4) 
£/ ~ ^(ri2 + r 34 ) + ^(ri 3 + r lA + r 23 + r 24 ) + J(rn + r i2 + r i3 + r i4 ) 

in ($(55)3(99)351^/1^(99)3(55)95)- ( 5 ) 

M° depends on the flux-tube topology and is denoted here as M^, M° 3 , M^--] and 
M[ qq ][ qq ]q f° r the configurations shown in Fig. 1(a), (b), (c), and (e), respectively, ([qq] and 
[qq] indicate (qq)z and (99)3, respectively.) 

In the present calculation, we ignore other terms such as tensor and spin-orbit interactions 
in the OGE potential, and we do not introduce flavor-exchange interactions. As shown later, 
the major properties of the normal hadron mass spectra is qualitatively reproduced by the 
present Hamiltonian. 



III. MODEL WAVE FUNCTION AND PARAMETERS 

We solve the eigen states of the Hamiltonian with a variational method in the AMD 
model space proposed in the previous paper. The AMD wave function in a quark model is 
given as follows. 



*(Z) = (1±P)A <t> Zl <t>z 2 ---<t>z Na ^ x 



( 1 ^ 3/4 

ex p 



V2bZ t f + \Z\ 



(6) 
(7) 



where 1 ± P is the parity projection operator, A is the anti-symmetrization operator, and 
the spatial part 4>z l of the ith single-particle wave function given by a Gaussian whose center 
is located at Zi in phase space. The spin function, $ 5 , is given as 

® S = E c mi ... mNq \mi- ■ -mjvjs, (8) 

mi,--,m Nq 

where \m) s (m =t, I) is the intrinsic-spin function. $ x is the flavor-color function. For 
example, the flavor-color function for the tetraquark udss system with color-configuration 
(99)3(99)3 i s written as 

$ x = \udss) <g> e abc e efc \abef) c - (9) 

In the present wave function we do not explicitly perform isospin projection, but the wave 
functions obtained by energy variation are found to be approximately isospin-eigen states 
in most cases due to the color-spin symmetry. 

As already mentioned, different flux-tube topologies appear in each of the q 2 q 2 and the 
q A q systems. Since the transitions between the different string configurations are of higher 
order in the strong coupling expansion, we ignore the coupling and perform a variational 
calculation within a single flux-tube topology. In the present calculations, we adopt only the 
connected flux-tube configurations given in Figs. 1(c) and (e), because we are interested in 
the confined and narrow states. This is regarded as a kind of bound state approximations. 

In the numerical calculation, the linear and Coulomb potentials are approximated by 
seven-range Gaussians. We use the same parameters as those adopted in Ref. [26]: 

a c = 1.05, 
A = 0.13 fm, 
m q = 0.313 GeV, 
m s = 0.513 GeV, 

a = 0.853 GeV/fm. (10) 
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(a) (b) 




FIG. 1. Flux-tube configurations for qq meson (a), q 3 baryon (b), q 2 q 2 states (c),(d), and q 4 q 
states (e),(f). For the q 2 q 2 states, the exotic tube configuration (c) corresponds to the (99)3(99)3, 
and the disconnected tube (d) represents the meson-meson state, (99)1(99)1. The configurations, 
(qq)^(qq)sQ and (999)1(99)1 for the q A q system are illustrated in figures, (e) and (f), respectively. 

Here, the quark-gluon coupling constant (a c ) and the string tension (a) are chosen so as to 
fit the mass splitting among N, A and AT* (1520). The width parameter (b) is chosen to be 
0.5 fm. 

IV. RESULTS AND DISCUSSIONS 

A. mesons and baryons 

In Fig. 2, we display the calculated masses of the conventional mesons and baryons 
compared with the experimental data. The zero-point energy of the string potential for the 
3q system is chosen to be M° 3 = 972 MeV to fit the nucleon mass, while for the qq is 
adjusted to be 584 MeV to reproduce the p-meson mass. It is shown that the systematics 
of the mass spectra are reasonably reproduced by the present calculations, except for the 
pseudoscalar mesons. 
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FIG. 2. Mass spectra for qq mesons and 3q baryons. The experimental and calculated masses 
are shown by dashed and solid lines, respectively. We adjust the zero-point energy of the string 
potential for the qq system as M qq = 584 MeV to fit the experimental /3-meson mass. For the 3g 
system, M° 3 is chosen to be M° 3 = 972 MeV to reproduce the nucleon masses. 

B. tf+-meson(/ = 0, J p = 1+) 

As mentioned before, the zero-point energy (M°) of the string potential depends on the 
flux-tube topology. We, here, phenomenologically deduce the unknown M[ 9? ][^] for the 4g 
system with the help of the systematics of M q g, M q z and M[ qq ^ qq jg for the normal meson, 
baryon and pentaquark systems. In a previous paper, we applied the present method to the 
uudds system and studied the properties of + . In the results, it was predicted that the 
three narrow states, I = 0, J p = l/2 + ,3/2 + states and the I = 1, J p = 3/2~ state may 
degenerate in almost the same mass region. M[ qq ][ qq ] q = 2375 MeV was chosen to fit the 
theoretical mass to the observed @ + mass. 

In Fig. 3, the adopted M qq , M q z and M^ qq ^ qq -^ q are shown as a function of the 
quark number (N q ). If the string potential is assumed to be a two-body linear-potential, 
— a s FiFj(r — r ), the potential is equivalent to the approximated flux-tube potential in 
Eqs. 2-5 with the relations a = |a s and M° = ^Nq. As a result, M° should be proportional 
to N q in the pairwise confining potential, which leads to the relation M^ qq ^ qq -\ q /M q 3 = 5/3. 
However, as shown in Fig. 3, the ratio M^ qq ^ qq -\ q /M q 3 is larger than 5/3 in the present 
model, and also in the pentaquark study with a constituent quark-model calculation in Ref. 
[27]. This means that we need an extra attraction in the [99] [99] 9 system in addition to 
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FIG. 3. Adopted zero-point energy (M°) of the string potential as a function of the quark 
number N q . M qq , M° 3 and MP are adjusted to reproduce the />meson, nucleon and 

+ masses, respectively. The M? qq ^--^ for the (99)3(99)3 is deduced by assuming the linear 
function M°(N q ) = aoN q + bo, where the parameters ao and 60 ar e determined by fitting 
(i) M°(N q = 2) = M$ q = 584 MeV and M°(N q = 5) = M [qq][qq]q = 2375 MeV, or (ii) 
M°(N q = 3) = Mj 3 = 972 MeV and M°(N q = 5) = M [qq][qq]q = 2375 MeV. The circle and triangle 
indicate the M? qq ^ qq ^ values obtained by the former fitting (i) and the latter one (ii), respectively. 

the pair- wise confining potential to understand the absolute mass of the + (1.54) within 
constituent quark models. We consider the dependence of M° on the tube topology as the 
many-body potential, and we here accept the value M[ qq ^ qq ] q = 2375 MeV [26] adjusted to 
the experimental 6 + mass. Then we phenomenologically determine the for 4g states 

by using the M?,,^ as an input as follows. 

We assume the linear function M°(N q ) = aoN q + bo and determine the parameters (ao, bo) 
by fitting the M values for the qq and [qq] [qq]q systems as (i) M°(N q = 2) = M qq - = 584 MeV 
and M°(N q = 5) = M[ qq ][ qq ] q = 2375 MeV. We also use another parameter set for (a , b ) by 
fitting the Mo values for the q 3 and [ag][gg](Z systems as (ii) M°(N q = 3) = M q3 = 972 MeV 
and M°(N q = 5) = M[ qq ][ qq ] q = 2375 MeV. With the obtained parameter sets a and b , we 
obtain M [qq]m = 1785 MeV and M [qq]m = 1679 MeV from M°(N q ) = a N q + b , (N q = 4) 
for the former fitting (i) and the latter one (ii), respectively. 

Now, we apply the flux-tube quark model with AMD to the (ud)s(ss) 3 system and 
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calculate the $ + (J p = 1 + ) mass. We use the above-determined zero-point energies, 
(i)M [ ° g][ -- ] =1785 MeV and (ii)Mg (?][ -- ] = 1679 MeV. The calculated #+(/ = 0,J P = 1+) mass 
is 1.37 GeV in case (i) and 1.46 GeV in case (ii). The results indicate that the $ + (J p = 1 + )- 
meson may exist around 1.4 GeV, near the KK* threshold (Fig. 4). 

Although we do not put a priori assumptions for spin and spatial configurations of 4 
particles, the -# + (J p = 1 + ) wave function obtained by the energy variation is dominated by 
the component with the spin-zero ( ud)% and the spin-one (ss)s in the spatially symmetric 
orbit, (Os) 4 . The spin-zero (ud)s gain the color- magnetic interaction, while only the spin-one 
configuration is allowed in the spatially symmetric (ss)z pair. Therefore this is consistent 
with the naive expectation in the diquark picture. 

We comment on the accuracy of approximations Eqs. 3,4 and 5 for the tube length in 
the Hamiltonian. As discussed in a previous paper [26], the tube length (Lf) is reasonably 
simulated by the approximated tube length L app given by Eqs. 3,4 and 5, while L app is exactly 
equal the Lf in the qq system. If we assume the harmonic oscillator (0s) Nq configurations and 
ignore the antisymmetrization of quarks, we can calculate the ratio of the L app to the exact 
tube length (L f ), which is denoted by L^/L ( f 0s) . The ratio L^/Lf s) is 0.91, 0.86 and 0.84 
for the g 3 , (qq) 3(^)3 and (qq)^(qq)^q systems, respectively. In order to examine the effect 
of this factor on the tetraquark mass, we scale the tube length as L/ ~ / L^J x L app 
and estimate the expectation value of the Hamiltonian for the present wave functions. With 
the obtained energies, we retune the M by fitting the p-meson, nucleon and pentaquark 
masses, and reexamine the $ mass. We then find that the modification of the tetraquark 
mass by the scaled Lf is slight; the -# + (J p = 1 + ) mass can decrease by 10 MeV for case (i) 
and 4 Mev for case (ii). 

C. f\- meson 

Here, we discuss the possibility of a f\(J p C = l ++ )-meson with the Aq component. Since 
its dominant decay mode should be KK*, it may have some analogy with the tetraquark 
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lines) shown in the left and right panels were obtained by using (i) M? 



[<?<?] [<n] 



1785 MeV and 



(ii) M^j--] = 1679 MeV, respectively. The mass of the f\ with the 4q component was obtained 
by calculations of the usus with J PC = 1 ++ . The dashed lines are experimental masses of the 
/i-mesons. The experimental KK* threshold energy is also displayed by dotted lines. 

$(J P C = 1 ++ ). If we ignore the qq annihilation, we can calculate the mass of the J PC =\ ++ 
(us) 3 (us) 3 state within the present framework in the same way as for the tetraquark ^-meson. 
The (us)s(us)3(J PC =l ++ ) state corresponds to an isoscalar /i-meson and an isovector a\- 
meson, which degenerate in the present Hamiltonian. We concentrate on the /i-meson in 
the present paper. By using the same zero-point energies (i) M^--] = 1785 MeV and (ii) 
M® qq ^ = 1679 MeV, the (^5)3(^5)3 state with J PC =1 ++ for the /i-meson is calculated to 
be 1.45 GeV in case(i) and 1.54 GeV in case(ii) (Fig. 4). 

There are various theoretical interpretations of scalar and axial-vector mesons as P-wave 
qq states, 4g states and hybrid qqg states. In the mass region 1~1 .6 GeV, three /i-mesons, 
/i(1285), /i(1420), and /i(1510) are known, though the /i(1510) is not well established 
[42]. In the P-wave qq state, two /x-mesons are expected to appear in this energy region as 
partners in the qq nonet. It is considered that the lower one is dominated by the light-quark 
component (nn = uu + dd), and the major component of the higher one is the ss state. In 
the standard interpretation, the lowest /i(1285) is regarded as the nn state. On the other 
hand, /i(1420) and /i(1510) are candidates for the partner of the /i(1285) in the qq nonet, 
but the assignment is not yet confirmed. 

In the constituent quark model calculation of qq systems [43], the masses of two 1 ++ 



12 



states in the P-wave qq nonet are 1.24 and 1.48 GeV. The theoretical mass spectra of the 
1 ++ qq states seems to be consistent with the experimental ones if /i(1510) is assigned to 
be a partner of the /i(1285) in the flavor nonet. This is consistent with the assignment in 
Ref. [44]. On the other hand, an alternative interpretation that the /i(1285) and /i(1420) 
are qq partners is claimed in Refs. [12,42,45]. 

These interpretations lead to an indication that one of /i(1420) and /i(1510) may be 
a non-gg meson, while the other can be understood as being partners of the /i(1285) in 
the conventional P-wave qq states. In the present calculation of the (us)${us)z state with 
J PC =\ ++ , the theoretical mass in case (i) 1.45 GeV seems to be consistent with the /i(1420), 
while the mass in case (ii) 1.54 GeV energetically corresponds to the /i(1510), as shown in 
Fig. 4. The present results suggest the assignment of a /i-meson in the 1.4~1.6 GeV mass 
region as the 4q state. 

D. Width of i9 + -meson 

As mentioned above, we suggest that the $ + (J p = l + )-meson may appear in the energy 
region ~1.4 GeV near the KK* threshold. The expected decay modes are KK* and KKw. 
The width for the KKtc decay can be enhanced by the broad resonance, k(800) state, in 
the scalar Kir channel. On the other hand, the phase space for the direct three-body decay 
is generally suppressed. In fact, the typical width of three-body decays is of the order 
of a several MeV in the cu-meson width, for example. Therefore, if the branchings into 
two-hadron decays, KK* and Kk(800), are small enough, the width should be narrow. In 
order to discuss the stability of the i? + -meson, we, here, consider only the two-hadron decay 
modes. We should note that the decay mechanism of the f) + (J p = 1 + ) may be analogous 
with that of the /i-meson with the usus(J PC = 1 ++ ) state, where the KK* and KKtt(Kk) 
are expected to be dominant decay modes. In fact, in the decay of /i(1420) and /i(1510), 
which are the candidates of the usus(J PC = 1 ++ ), as mentioned before, the KK* and/or 
KKn modes are experimentally seen, and the former is dominant in the /i(1420) decay 
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modes [42]. For decay into the KK*, the S'-wave channel is allowed, while the Kk decays 
should be a P-wave. 

First, we give a rough estimation of the $ + width for the KK* decay by assuming that 
the coupling for fi — > KK* (g^xR*) is the same as that for i9 + — > KK* (g^KK*)- The width 
is approximated by the product of the coupling and the phase space. We take into account 
only the S'-wave decay, and estimate the phase space for the KK* decay by the imaginary 
part of the one-loop self-energy integral I(p) for the scalar mesons, 



Im[7(p)] = ^2 Im 



—ai ln(mj) — a 2 m(m 2 J — q In 



(11) 



(q - ai)(g - a 2 ) 
where, 

_ 1 /" i _ (mi+m 2 ) 2 \ / -. _ (mi-m 2 ) 2 \ 2 

q — 2 l v p 2 / v p 2 A 
fll = i(i-=^) 

a 2 = Kl-^)> (12) 

and mi and m 2 are the masses of the daughter particles. We take the rest frame p = (p , 0). 
As is well known, if the mi and m 2 are real values, Im[7(p)] is equal to the usual phase 
space q above the threshold (po > m i + m 2 ), and is zero below the threshold (|m 2 — mi| < 
Po < m i + WI2). However, since the daughter particle K* has a width of T^* ~ 50 MeV for 
Kn decays, the phase space should be finite, even below the threshold. In order to estimate 
the effect of V K * on the phase space, we take the masses of the daughter particles to be 
mi = fflp — iT 'k*/2 and m 2 = m^ in Eq. 12, as is done in Ref. [46]. We use the masses of 
K and K* as m K = 0.495 GeV and m K * = 0.892 GeV, and evaluate the width of the parent 
particle as Y = 2g 2 Im[/(m )], where mo is the parent particle mass and g is the coupling. 

We show the mass m dependence of the width r(m ) for the KK* decays in Fig. 5, 
where the coupling g is adjusted to reproduce the experimental full width Fj 1 ( 1420 ) = 56 
MeV at m = 1.42 GeV, and is assumed to be energy independent. If we adopt the case (i) 
calculation and the assignment of the /i(1420) as the 4q state, the width of the -$ + (l + ) is 
estimated to be T$ ~ 20 MeV. When we assign the /i(1510) as the Aq state and choose 
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the coupling constant g to fit the width 100 MeV of the /i(1510) at m = 1.51 GeV, we 
obtain ~ 80 MeV based on the case (ii) calculation. Here, we adopt the upper limit of 
the /i(1510) width in Ref. [42], though there still remains an ambiguity in the experimental 
data. 

Next, we discuss the effect of Kk decay on the width. In the above estimation of r#, 
the effect of this mode is already included in the total width of the /i-meson. In contrast to 
the S'-wave decay in the KK* channel, the 1 + states decay into the P-wave Kk state. Even 
though the P-wave decay is unfavored compared to the S'-wave decay, we should consider 
its effect, because the Kk threshold (~ 1.3 GeV) is lower than the KK* threshold (1.39 
GeV). We again assume that the coupling for fi — > Kk is the same as that for $ + — * Kk, 
and estimate the width by using the ratio of the phase space for the tf(m = 1.37 GeV), q#, 
to that for the /i(1420)(m = 1.42 GeV), q^, as T# = Tf x x q^/q^. Since k(800) is a broad 
resonance with V = 300 ~ 800 MeV, we may not apply the previous method by the single- 
pole approximation of the propagator in the one-loop diagram to estimate the phase space. 
Instead, we consider the phase space for the P-wave decay, \m — m K —m K \ 3 / 2 , and take into 
account the broad mass range of m K , where the m K is the k mass and is a real value. When 
we adopt the case (i) calculation, we obtain the ratio (qv/q^) of the phase space for the 
$(m = 1.37 GeV) to that for the /i(1420)(m = 1.42 GeV) as 0.5 for m K = 800 MeV. The 
ratio qo/qfi is 0.7 for the lower limit m K = 640 (Kn threshold) MeV, while q^/q^ = for 
m K > 900 MeV. On the avarage, the ratio q^/q^ in the Kk decay is almost the same as that 
in the KK* decay. Even if the branching ratio of the Kk(Kk) mode is 100% in -$ + (l + )(/i), 
an upper limit of T$ = 40 MeV is obtained from the = 56 MeV and the lower limit 
of m K . Also, in calculation (ii), the phase space ratio {q-a/q^) for /i (1510) — > Kk to that 
for the d — > Kk in the mass range m K > 640 MeV is calculated as q^/q^ < 0.8, which is 
consistent with that in the KK* decay. As a result, it is concluded that the ^-meson width 
is estimated to be O(20 — SOMeV"), no matter how great are the branching ratio in the KK* 
channel and that in the Kk channel. 

In the above discussion, we roughly evaluate the width of the -# + (l + )-meson by assuming 
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k (GeV) 

FIG. 5. Energy dependence of the two-meson KK* (892) decay width V. The coupling is chosen 
to fit the full-width T = 55 MeV of the /i(1420) at mo = 1.42 GeV, and is assumed to be energy 
independent. 

the same coupling for -# + (l + ) — > KK* and -# + (l + ) — > Kk as those for f\ — > fT^* and 
/i — > fT/t, respectively. In the present estimation we choose the couplings to fit the full 
width of f\. In the case that the couplings are enhanced via the annihilation and creation of 
a qq pair in the /i-meson, the width of the -$ + (l + )-meson might be smaller than the present 
estimation. 

V. SUMMARY AND DISCUSSION 

We discussed the possibility of the J p = 1 + state of the isoscalar tetraquark (S=+2), -# + - 
meson, with the udss content. If the pentaquark 6 + (1540) has the (qq)i(qq)?,q configuration, 
the $ + (J p = 1 + ) is expected to exist at an energy lower than, or close to, the pentaquark 
+ (154O) mass. This leads to a possible appearance of the -# + -meson as a resonant state, 
which may be observed in the K + K + n~ channel. 

We investigated the f) + (J p = 1 + ) with a constituent quark model. The flux-tube quark 
model with antisymmetrized molecular dynamics(AMD) was applied to 4q systems in the 
same way as in the pentaquark study in Ref. [26]. Based on the picture of a flux-tube model, 
we solved the 4q dynamics in the (qq)i(qq)3 model space by the variational method. The 
results suggest that the i) + (J p = 1 + ) may exist at around 1.4 GeV. Since the predicted 
mass of the d + (J p = 1 + ) is close to the lowest (KK*) threshold in the allowed two-hadron 



16 



decays, the present results imply that the -# + (J p = 1 + ) may exist as a resonance, and its 
width may be not very broad. 

We also calculated the J PC = 1 ++ (us) 3(^5)3 state, which is associated with a non-gg fi- 
meson. The calculated mass of the f\ with the 4g configuration suggests an interpretation 
that the Aq state may correspond to one of the /i-mesons in the 1.4 ~ 1.6 mass region. 
/i(1420) and /i(1510) are candidates of the 4g state. If we assume that the couplings 
for /1 — > ifF and /1 — > X/c* are the same as those for § + (J p = 1 + ) — > XX* and 
•d + (J p = 1 + ) — > i^K, respectively, we can evaluate the width of the f> + to be 20-80 MeV 
from the phase space for these two-hadron decay modes. Provided that the coupling for 
the two-hadron decays in the d + (J p = 1 + ) is small enough, the dominant decay should 
be a direct three-hadron decay, $ + — > KKir, with a small phase space, and hence the 
$ + (J p = 1 + ) should have a small width. 

Recently, the $ + -meson was discussed by Burns et al. [34] and by Karliner and Lipkin 
[47]. In Ref. [47], it was mentioned that the $(J P = + ) can be narrow, because the lowest 
allowed decay mode is a four-body KKnn channel with a small phase space. However, 
the J p = + state is forbidden in the isoscalar udss system within spatially symmetric 
configurations, and hence, the ■&( J p = + ) is expected to be energetically unfavored. Burns 
et al. predicted the d + (J p = 1~) with L = 1 at ~ 1.6 GeV with a width of O(10 - 100) 
MeV, which can decay into K + K°. Although the color-magnetic attraction may be larger 
in the $(J P = 1~) state than in the d(J p = 1 + ) state, the $(J P = 1~) must have the 
L — 1 excitation energy. Another claim for the d + (J p = l - ) state is that it can decay into 
P-wave KK states. The centrifugal barrier may not be high enough to stabilize the state 
much above the threshold energy. Therefore, we consider that the -d + (J p = 1 + ) is a better 
candidate of narrow tetraquarks. 

Our calculations of the tetraquark d + (J p = 1 + ) and the pentaquark 6 + are based on the 
color-configurations (qq)^(qq)^ and (qq)3(qq)?,q- We ignore (qq)e configurations because the 
(qq)e is an excited configuration, and is unfavored in the strong-coupling picture. We should 
remark that the preference of the spin-parity J p = 1 + in the $ + system does not depend 
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on the color configuration. This is because, in the isoscalar udss with a spatially symmetric 
configuration, other spin parities are forbidden, and therefore the spin parity is uniquely 
determined to be J p = 1 + . This means that, even in the color-configuration (ud) 6 (ss)g, 
as suggested in Refs. [34,47], only J p = 1 + is allowed in the spatially symmetric d state. 
If the 6 + has a triquark-diquark structure ((ud) Qs) 3 (ud)^ with the relative P-wave motion 
as proposed by Karliner and Lipkin [24], it is expected that the r d{J p = 1 + ) mass with a 
spatially symmetric (-uc?) 6 (ss)6 configuration would be close to, or could be smaller than, 
the + mass due to the kinematic energy gain. Another important channel is the meson- 
meson {qq)i{qq)\ configuration, though it can be expressed by a linear combination of the 
(ud)^(ss) 3 and (ud) 6 (ss)Q configurations. If the f> + (J p = 1 + ) mass is smaller than the lowest 
meson-meson threshold KK* , mixing of the (<?<?) 1 (<?<?) 1 state may not have a significant effect 
on the $ width, because the q 9 + {J p = 1 + ) is bound in the KK* channel. In the case that 
the -# + (J p = 1 + ) is heavier than the KK* threshold, channel mixing should be taken into 
consideration to estimate its width. In IV D, we gave an analysis of the widths by assuming 
that the decay mechanism and the coupling in the i) + (J p = 1 + ) are the same as those in the 
/x-meson. Within this assumption, the channel-coupling effect on the widths is effectively 
included in the adopted total width of the /i-meson. For a further detailed discussion of 
the stability of the -#(J P = 1 + ), coupled calculations of different color configurations are 
required. 

We should point out that the allowed decay channels are different among these three 
predictions (J p = 1 + , 1~, and + states) of the -$ + -mesons. For the -$ + (J p = 1 + ), d + (J p = 
and -# + (J p = + ), the decay modes are KKn, KK and KKirn, respectively. In order 
to establish the udss content in the tetraquark -# + -meson, it is necessary to observe not K°, 
but at least two K + s, because the K° contains the sd component as well as the ds. In that 
sense, the K + K + n~ decay from the i?(J p = 1 + ) predicted in the present work is suitable 
for an experimental tetraquark search. 

In the constituent quark model calculations of Refs. [4,6], there is no indication for the 
existence of multiquark hadrons, except for the /o(980) and a (980) as KK molecules [4], 
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while interpretations of scalar mesons, like / (600), k(800) and D s j(2317), by four-quark 
states have been suggested in several quark model calculations [1-3,48-51]. In Ref. [6], the 
absence of the pentaquark was claimed, which is inconsistent with the observation of the 
pentaquark 6 + . In order to explain the absolute mass of the 6 + (1540) within constituent 
quark models, we need an extra attraction for the multiquark system in addition to a pair- 
wise interaction. In the present calculation, we did not obtain a quantitative value of this 
extra attraction from fundamental theory; instead, the many-body potential was taken into 
account in terms of the flux-tube potential, and we phenomenologically evaluated it by using 
the observed 6 + mass as a input. 

In the present work, we used a simple Hamiltonian with the confining force, and the 
Coulomb and color-magnetic terms in the OGE potential. For a systematic description 
of the hadron spectra, there still remain such problems as fine tuning of the interaction 
parameters and inclusion of the tensor and spin-orbit interactions in the Hamiltonian. To 
obtain theoretical insights into multiquark hadron physics, an experimental search for the 
tetraquark is necessary as well as further experimental studies on the pentaquark. We 
conclude that the d + (J p = l + )-meson is proposed as a good candidate of the tetraquark, 
which would be observed in the K + K + tt~ decay channel. 
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